Dynamical and bursty interactions in social networks 
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We present a modeling framework for dynamical and bursty contact networks made of agents in 
social interaction. We consider agents' behavior at short time scales, in which the contact network is 
formed by disconnected cliques of different sizes. At each time a random agent can make a transition 
from being isolated to being part of a group, or viceversa. Different distributions of contact times 
and inter-contact times between individuals are obtained by considering transition probabilities with 
memory effects, i.e. the transition probabilities for each agent depend both on its state (isolated 
or interacting) and on the time elapsed since the last change of state. The model lends itself to 
analytical and numerical investigations. The modeling framework can be easily extended, and paves 
the way for systematic investigations of dynamical processes occurring on rapidly evolving dynamical 
networks, such as the propagation of an information, or spreading of diseases. 

PACS numbers: 89.75.-k,64.60.aq,89.65.-s 



Recently, technological advances have made possible 
the measure of social interactions in groups of individu- 
als, at several temporal and spatial scales and resolutions 
showing that human activity obeys scaling law and sta- 
tistical features which reveal long time correlations and 
memory effects. Evidence comes from data on email ex- 
changes [H-Q , mobile phone communications 0, IH , spa- 
tial proximity 6--9L web browsing [10( and even face to 
face interaction [Ill4l3|. In this respect the traditional 
framework of models used for risk assessment and com- 
munication, which describe human actvity as a series of 
Poisson distributed processes, need to be changed in fa- 
vor of new models which take into account the occurrence 
of burstiness in many aspects of human activity. 



Social interactions give rise to social [14J, [iff and col- 
laborative jig networks characterized by a complex evo- 
lution. In these networks, links are constantly created 
or terminated, and the social network of an individual 
evolves at different levels of organization. After the pio- 
neering papers on complex networks showing that many 
social networks are small world and display heteroge- 
neous degree distributions [17| . and that these network 
topologies strongly influence the dynamics taking place 
on the networks |18| . a number of papers have been de- 
voted to modeling the dynamics of social interactions. 
Issues investigated in this context are in particular com- 
munity formation and the evolution of adaptive 
dynamics of opinions and social ties 2l|- 24 1 . 

The evidence coming from the analysis of social contact 
data calls for new frameworks that integrate these models 
with the bursty character of social interactions. The du- 
ration of contacts between individuals or groups of indi- 
viduals display indeed broad distributions, as well as the 
time intervals between successive contacts @, [H, 12, 2o| . 
Such heterogeneous behaviors have strong consequences 
on dynamical processes [H, [2(| , and should therefore be 



correctly taken into account. It is therefore necessary 
to introduce this fundamental aspect on human activ- 
ity in models of social interactions, possibly reconstruct- 
ing then social networks by aggregating the network of 
contacts over a certain period [26l - [28| . The mode ling 
literature in this area being still in its infancy (2|| I29I - 
[32)], it is important to develop simple, generic, easily im- 
plementable models of dynamical networks which repro- 
duce the empirical facts observed in contact duration and 
inter-contact intervals. 

In this Letter, we take a step in this direction, focus- 
ing on short timescales such as the ones involved when 
people interact in social gatherings (e.g., scientific confer- 
ences). We define a simple agent-based model for rapidly 
evolving sparse dynamical networks, aimed at describing 
the dynamics of human social interactions in the con- 
text of small discussion groups. In particular we are in- 
terested in investigating basic mechanisms which could 
be responsible for various contact duration distributions. 
The model is kept simple, so that it can be easily sim- 
ulated. It is accessible to analytical investigations in a 
certain number of cases. It can also be easily extended 
or modified. For instance, the population of agents is 
considered homogeneous (i.e. every agent is assumed to 
have the same dynamical parameters) and an extension 
to heterogeneous populations can easily be envisioned. 

The dynamical network under study is formed by dis- 
connected groups of agents which evolve by successive 
mergings and splittings. In particular at each time step 
an agent can either leave or remain in its group, or intro- 
duce an isolated agent to its group. The general formula- 
tion of the model allows to describe a variety of behaviors 
of the dynamical networks. In particular, the duration 
of contacts between individuals can display either nar- 
row or broad distributions. A narrow distribution is for 
instance obtained by simply assuming that each agent 
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leaves a group or invites a new agent in its group with 
a time-independent probability. On the contrary, broad 
distributions of contact durations, 



similar to those ob- 
served in empirical studies @-H, 11, [25|, are obtained 
through a reinforcement dynamics of the interaction, that 
can be summarized as "the longer an agent interacts with 
a group, the less it is likely to leave the group; the more 
the agent is isolated the less likely it is to interact with 
a group". This dynamics, reminiscent of the preferen- 
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tial attachment in the context of complex networks [33J , 
could be argued to stem from Hebbian-like mechanisms 
at the underlying cognitive level. In general, for both nar- 
row and broad distributions of interaction times, larger 
groups are found to be less stable than smaller ones. This 
is also observed in the data [ll| and can be simply ex- 
plained: the lifetime of a group depends on the decisions 
of all its members. In a first approximation these deci- 
sions correspond to independent events, therefore groups 
with more agents become less stable. Interestingly, our 
model also exhibits a dynamical transition towards the 
formation of large size group. This transition, supported 
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by some measurements in animal behavior 
not observed in human behavior, and corresponds thus 
to parameter values where the model loses its applicabil- 
ity to the description of human social interactions. Note, 
nevertheless, that the formation of large social organiza- 
tions and cities, demonstrate that in humans, large group 
formation occurs at a different level of organization. 

The model we propose considers a fixed population of 
N agents, interacting in a limited space, as for example 
in a conference venue Therefore, in a first 

approximation we neglect the spatial dispersion of the 
agents and assume a well mixing dynamics. Each agent 
can either be isolated or belong to a group with other 
agents, and the groups define an instantaneous contact 
network. During the dynamics, agents can join other 
agents or on the contrary leave the group they belong to. 
More precisely, each agent i is characterized by two vari- 
ables: the number pi of other agents with which it is in 
contact (i.e. its degree in the network) and the time ti at 
which pi last evolved. At each time step t, an agent i is 
chosen at random. If i is isolated (pi = 0), i changes its 
state with probability bof{t, U). In this case, another iso- 
lated agent j is chosen with a certain probability H(t, tj), 
and i and j form a pair (pi —> 1, pj —> 1 and U — > t, 
tj — > t). If on the other hand i is part of a group Q of 
size greater than one (i.e. i has degree p > 0), a change 
of state occurs with probability b\f(t,ti). When this oc- 
curs, agent i can either leave the group (probability A) 
or introduce an isolated agent in the group (probabil- 
ity 1 — A). If i leaves the group Q and becomes isolated, 
Pi — > , and pj — >• p — 1 Vj 6 Q \ { i } , and as a consequence 
of this event the time of the modified nodes is reset to 
t, i.e. tt — > tW € Q. Hi introduces to the group an 
isolated agent j, chosen again with probability H(t,tj), 
then pi — >• p+lVl 6 QU {j} and each agent £ in Q U {j}, 



changing state at time t, sets ti — > t. The parameters 
bo and b\ determine the tendency of the agents, respec- 
tively isolated or in a group, to change their state, while 
A controls the tendency either to leave groups or on the 
contrary to make them grow. The model's dynamical 
behavior depends also on the functions / and II. 

In order to make contact with empirical data, the main 
quantities of interest concern the time spent by agents in 
each state, the duration of contacts between two agents, 
and the time intervals between successive contacts of 
an agent. We can gain insight into these properties by 
writing rate equations for the evolution of the number 
N p (t, to) of agents which are at time t in state p since to'. 
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where in the sums t' < t. These equations can be sim- 
plified and solved in certain cases, and the distribution 
Pp(r) of (normalized) times t — (t — to)/N in which an 
agent remains in a given state p can then be deduced. 
Let us for definiteness assume that / and II are station- 
ary functions so they depend only on t — to', it is then 
natural to look for a stationary state, reached at large 
enough times, such that a, r, (7r Pjg ) are constants, and 
N p (t,to) = N p (t — to). If for instance / is a constant, it 
is easy to see that the {N p } eN decay exponentially with 
time, so that the {Pp( T )} p> o are as weu exponentially 
decaying functions. 

We consider the more interesting case of / and II de- 
caying with t — to', the more an agent is in a state, 
the less probable it becomes to change state, as previ- 
ously described in the self-reinforcement mechanism. For 
sake of simplicity, we focus on the case / = II, so that 
r = 1, which allows to simplify the computations [36[. 
Computations can be carried out completely for instance 
in the case f(t — t ) = H(t — t ) = (1 + t) _1 , where 
t = (t — to)/N. The choice of this scaling is consis- 
tent with the scaling of email communications and other 
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FIG. 1: Average state (p) in the stationary state for different 
sets of parameters and a population of N = 2000 agents, (a) 
bi = A = 0.7; (b) & = A = 0.7; (c) b = 6i = 0.7. In (c) we 
observe the divergence of the average group size as A — > 0.5 + . 
The lines show Eq. ©. 



human activity @, E3]- In particular, N p (t) for p > 1 
is readily seen to decay as a power-law with exponent 
(p + More involved computations are needed to ob- 
tain the decay exponent of Nq. Writing {7r Pi9 } as func- 
tions of {N p } GN , / and II, we can obtain recurrence 
relations for (to) p>0 an d deduce a — (2A — 1)~ 1 , so 
that 

P (r) oc (l + T)- 1+b «^ (3) 
P p (r) oc (1 + t )-i+(p+i)6i , P >1 . (4) 

The previous analytical results are obtained under the 
conditions bi > 1/2, A > 1/2, b > (2A - 1)/(3A - 1), 
which determine the phase diagram of the model: out- 
side these boundaries, the hypothesis of stationarity is 
violated. 

It is also possible to compute the average state of the 
agents in the stationary state, as 
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where 



To = 
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For A — > 0.5 + the average group size (p) diverges indicat- 
ing that, in this limit, the non-stationary state is domi- 
nated by the formation of a large group of size O(N). 

We have performed numerical simulations of dynami- 
cal networks generated by the present model, with dif- 
ferent /, II, values of the parameters bo, bi, A, and sizes 
N. We will here show the simulations corresponding to 
/ = II = (1 + r) _1 , in order to compare with the analyt- 
ical predictions presented above. We first show in Fig. Q] 
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FIG. 2: Distribution P P {r) of times during which an agent 
does not change connectivity p. N = 10000, &o = &i = 0.7, 
A = 0.8, and the simulation is run for T = 10 N time steps. 
The lines are the analytical predictions. 



the average agent state as a function of the different pa- 
rameters, recovering the behavior predicted in Eq.([S]). 
The average state increases with bo, decreases with A, 
and presents a non-monotonous behavior with b\ . Figure 
[H displays the distributions {-P p } < p <4 01 time spent in 
the various states. These distributions are power-laws, 
in perfect agreement with the analytics. We also note 
that, for / = II = (1 + T )-" with v ± 1, {P P (r)} p£N 
can be shown analytically to become either stretched ex- 
ponentials (y < 1) or power-laws (y > 1), and we have 
also checked this behavior numerically. The broadness of 
the distributions is therefore not limited to the particular 
case described above, but is quite robust with respect to 
changes in the microscopic rules. 

In Fig. [31 we also show the distribution of contact du- 
rations between two agents (which is different from Pi: 
two agents remain in contact when they are joined by a 
third, but leave the state p = 1), of triangle durations, 
and of the time intervals between the starting times of 
two successive contacts [ll|. This last quantity is highly 
relevant in the context of causal processes, as it gives the 
time scale on which an agent can propagate an informa- 
tion or a disease after receiving it. All these distributions 
are broad, similarly to empirical observations [ill |25| . 

Let us finally mention that, when considering parame- 
ter values outside the validity of the stationary state an- 
alytical computations, different scenarios are observed, 
depending on A: if A > 0.5, the average state slowly de- 
creases (towards if bo < 0.5, and 1 if b\ < 0.5) while, for 
A < 0.5, a large cluster appears, with size is proportional 
to N, and lasting on a diverging timescale. Interestingly, 
even in this non-stationary case, the shape of the distri- 
butions P p (r) may remain stationary (not shown). This 
is particularly relevant as most empirical data are neces- 
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FIG. 3: Distributions of (a) duration of a contact between two 
agents; (b) time intervals between the beginnings of successive 
contacts of an agent A with two different agents B and C; (c) 
duration of a triangle, bo = 0.7, A = 0.8, bi = 0.6 and 0.9. 
N = 1000, T = 10 5 N. 



sarily obtained in non-stationary environments. 

In this Letter, we have proposed a modeling frame- 
work for dynamical networks in the context of interact- 
ing social agents. Both broad and narrow distributions 
can be obtained, corresponding to different social situa- 
tions. The present framework can be developed in sev- 
eral research directions. First, many variations of the 
microscopic rules may be thought of and implemented, 
in order to model more precisely mechanisms of social 
contacts in various contexts or even of animal behavior. 
For instance, merging and splitting of groups could be in- 
troduced, as well as heterogeneity between agents to take 
into account different propensities to interact or to cre- 
ate groups. Moreover, it will be interesting to investigate 
how the properties of the interaction durations shape the 
resulting aggregated networks on various timescales. Fi- 
nally, model dynamical contact networks can be used as a 
support for the simulation of dynamical processes taking 
place on dynamical networks, such as information spread- 
ing in a conference: the spreading process, although tak- 
ing place on an extremely sparse network which is at any 
time formed of disconnected groups, may overall concern 
the whole population of agents, thanks to the dynam- 
ics of the agents who move from one group to another. 
The fact that the various network characteristics (such 
as the broadness of the distribution of contact durations 
and inter-contact times) can be controlled by changing 
the model's parameters will then make it possible to un- 
derstand better the effect of these characteristics on the 
dynamical processes under scrutiny. 
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